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Abstract. Recent discovery of D s states suggests existence of radial excitations. Our semirelativistic quark 
potential model succeeds in reproducing these states within one to two percent of accuracy compared with 
the experiments, D s o(2860) and (2715), which are identified as + and 1~ radial excitations (n = 2). We 
also present calculations of radial excitations for B/B a heavy mesons. Relation between our formulation 
and the modified Goldberger-Treiman relation is also described. 
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1 Introduction 

BaBar has recently announced the discovery of a new D s 
state, which seems to be cs state, pQ D s0 (2860). Subse- 
quent to this experiments Belle has observed a new state 
of -D*(2715) whose spin and parity is determined to be 

We were the first in predicting + and 1 + states of 
D s and D particles, D s0 (2317), A,i(2460), L»5(2308), and 
Z)^(2427),[3], and we have also succeeded in reproducing 
higher resonances of B/B s particles, Bi(5720), £?|(5745), 
and B* 2 (5839),[H using our semirelativistic quark poten- 
tial model. Our model succeeds in lowering + and 1 + 
states of D s so that these mass values are below DK/D*K 
threshold which other models do not easily succeed. This 
model respects both heavy quark symmetry and chiral 
symmetry in a certain limit of parameters, [5] which relates 
our formulation with the idea of the modified Goldberger- 
Treiman relation proposed in Refs. [HE]. Hence it is nat- 
ural to try to explain the newly discovered D s states by 
using our semirelativistic model and we are again success- 
ful in reproducing these states discovered by BaBar and 
Belle. 

To interpret the state D s o(2860), there are arguments 
that this cs state is explained to be a scalar by a coupled 
channel model[S], or that it is a J p = 3~ state[9], or that it 
can be explained by using a phenomenological interaction 
terms like our quark potential model, [10) or that it can be 
analyzed by using the 3 Po model. [TT] 

Starting from the astonishing discovery of D s j parti- 
cles with narrow decay width by BaBar and CLEO, and 
confirmed by Belle, a series of successive experiments on 
the spectrum of a heavy-light system, i.e., heavy mesons, 
heavy quarkonium, and heavy baryons, stimulates theo- 



rists to explain all these spectra as well as their decay 
modes. See the recent reviews of Refs. [j"2"llT5] . It seems 
that a new era of spectroscopy is opening, which is chal- 
lenging to theorists to solve these spectra at the same 
time. 



2 Numerical Calculation 

Our model starts from a Hamiltonian with a scalar confin- 
ing potential together with a Coulombic vector potential, 
we expand the whole system, i.e., Hamiltonian, wave func- 
tion, and eigenvalue is expanded in 1/toq, and we solve 
equations order by order consistently. In the actual nu- 
merical calculation, we expand a wave function in a power 
series of relative coordinate with some weighting exponen- 
tial times power functions. A wave function has positive 
components of a heavy quark due to the lowest order con- 
straint when expanding in 1 /toq and four components of 
a light antiquark. Thus the lowest eigenfunction is two by 
four matrix. Other than angular part, its radial part can 
be written like below. 



Ufc(r), v k (r) - w k (r) (^j 



exp 



-{m q + b)r--[- 



where A; is a quantum number of an operator, — (3 q (S q ■ L + 
with S q light quark spin and L light quark angular mo- 
mentum, which distinguishes uniquely each state [3J, for 
instance, k = —1 for a spin multiplet (J p = — 
k = +1 for (0 + ,l + ), etc., Ufc(r) and Ufc(r) are upper 
and lower components of radial wave functions, and A = 

*J k 2 — (4cv s /3) 2 , where a and b are included in a scalar 
potential as S(r) = r/a 2 + b, a s is a strong coupling, 
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Table 1. Most optimal values of parameters. 



Parameters 


„ n—'Z 


a (GeV -1 ) 


b (GeV) 




0.344 


1.939 


0.0749 



m U}d (GeV) 


m a (GeV) 


m c (GeV) 


m b (GeV) 


0.0112 


0.0929 


1.032 


4.639 



and Wk(r) is a finite series of a polynomial in r, Wk{r) = 
J^i^Q 1 a i ( r / a Y ' which takes different coefficients for Ufe(r) 
and Vk(r). In actual calculation, we have used N = 7. 
Hence we can in principle obtain seven different radial ex- 
citations. 

In this paper, we calculate the most optimal values of 
parameters so that recently discovered and known n = 
2 (the first radial excitation) D s particles are all fitted 
well around one percent of accuracy compared with the 
experiments. Here only a strong coupling a s is modified 
and other parameters are kept the same as in [4] , in which 
we have obtained a s = 0.261 both for D and D s . These 
are presented in Table Q] at the first order of calculation 
in p/rriQ with p being internal quark momentum and toq 
heavy quark mass. 

With these values of parameters, we obtain n = 2 
masses of D s and D at the same time. The results are 
shown in Table [3] for D s . We also predict n = 2 states 
for D, B, and B s states, which are shown in Tables HI [5j 
and [5] assuming the same strong coupling a s for D s , which 
may actually be different for B/B s particles. In these Ta- 
bles, pi and rii are i— th order positive and negative com- 
ponent contributions of a heavy quark, respectively, and 
Ci — Pi + Jij. When one carefully looks at these Tables, 
one notices that values of higher states 3 Di and " 3 D2 V 
are not reliable even though we have listed in Tables for 
consistency with the former calculations. 



3 Mass Gap 

Noticing the mass gaps between spin multiplets, (CP, 1~) 
and (0 + ,l + ) for D s j mesons, are almost equal to each 
other, the modified Goldberger-Treiman relation is pro- 
posed by Bardeen et al.[Hl[7] to understand the facts, i.e., 
the underlying physics might be chiral physics. In other 
words the mass gap is caused due to the chiral symme- 
try breakdown which is expressed by this relation. Fur- 
ther they have assumed hyperfine splittings due to break- 
down of heavy quark symmetry (inclusion of 1/mg cor- 
rections) are the same within two spin multiplets, (CP, 1~) 
and (0 + , 1 + ), so that the mass gap is not affected by this 
hyperfine splitting. 

In our formulation this is explained in Fig. [T] 5j Heavy 
quark symmetry reduces the original Hamiltonian into the 
one without spin structure after projecting wave functions 
into positive and negative components of a heavy quark. [3l 
[5] Chiral symmetry of a light quark is realized by taking 
a limit of m q — > and S(r) —* 0, in which case all the 
members of two spin multiplets, (CP,1~) and (0 + ,l + ), 



Table 2. Theoretical mass gap. Values in brackets are exper- 
iments. Units are in MeV. 

Mass gap (n = 1) D U s B B s 

0+ - 0" 414 (441) 358 (348) 322 239 

1+ - 1" 410 (419) 357 (348) 320 242 



[n 


= 2) 


D 


D 3 


B 


B s 


0+ 


-0" 


308 


274 


206 


160 


1+ 


- P 


350 


327 


216 


171 



are degenerate. When the light quark mass and a scalar 
potential are turned on, then the degeneracy due to chiral 
symmetry is resolved and the mass gap corresponding to 
the modified Goldberger-Treiman relation is given by 

AM = M {k = +1) - M (k = -1), 

where Mo(fc) is a degenerate mass for a quantum number 
k, which appears in Tables [3H6] to distinguish states. This 
is described as AM — g^f^/GA in 0. They have assumed 
dominant interaction terms for hyperfine splitting due to 
1/rriQ corrections so that mass gaps between + — CP and 
1 + — l - are almost equal to each other, which seems to 
hold in our formulation. Our hyperfine splittings due to 
l/m,Q to this mass gap AM are calculated and we have 
added those to degenerate mass gaps between + — CP 
and 1 + — 1~ for D, D s , B, and B s with n = 1, which 
are given in Table Our dynamical calculation supports 
the assumption that the mass gaps between + — CP and 
1 + — 1~ are almost equal to each other in the case of n = 1. 
In the second row of the same Table values for n = 2 are 
given, which is not as good as the case for n = 1. 
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Fig. 1. Procedure how the degeneracy is resolved in our model. A quantum number k is defined in the main text. 







Table 3. D s (n 




2) meson mass spectra (first order). Units 


are in MeV. 




• 2s+i Lj{J p ) 


Mo 


ci/Mo 




pi /Mo 


m/Mo 


M ca l c 


M obs 


L So{0~) 


2328 


1.006 xlO" 


i 


0.919 xlO" 1 


8.695 xlO- 3 


2563 




3 S 1 (1~) 




1.830 xlO" 


i 


1.824 xlO" 1 


5.744 xlO" 4 


2755 


2715 


3 Po(0+) 


2456 


1.553 xlO" 


i 


1.470 xlO" 1 


8.245 xlO" 3 


2837 


2856 


" 3 Pi"(l+) 




2.551 xlO" 


i 


2.543 xlO" 1 


7.667 xlO" 4 


3082 




" 1 Pi"(l+) 


2585 


1.969 xlO" 


i 


1.966 xlO" 1 


2.531 xlO" 4 


3094 




3 P 2 (2+) 




2.209 xlO" 


i 


2.209 xlO" 1 


5.070 xlO" 7 


3157 




3 Dr(l-) 


2391 


8.605 xlO" 


i 


8.600 xlO" 1 


5.016 xlO" 4 


4449 




" 3 P 2 "(2") 




-4.287 xlO" 


i 


-4.287 xlO -1 


5.482 xlO" 7 


1366 





Table 4. D(n — 2) meson mass spectra (first order). Units are in MeV. 



' 2s+i Lj(J p ) 


Mo 


ci/Mo 


Pi /M 


ni/Mo 


Mcalc 


M ob . s 


L S (Q-) 


2241 


1.078 xlO -1 


0.975 xlO" 1 


1.038 xW~' 2 


2483 




; %(1") 




1.917 xlO" 1 


1.910 xlO" 1 


6.882 xlO" 4 


2671 




3 P (0+) 


2418 


1.540 xlO" 1 


1.444 xlO" 1 


9.621 xlO" 3 


2791 




" 3 Pi"(l+) 




2.493 xlO" 1 


2.488 xlO" 1 


5.352 xlO" 4 


3021 




" 1 P 1 "(1+) 


2491 


2.076 xlO" 1 


2.070 xlO" 1 


5.956 xlO" 4 


3008 




3 P 2 (2+) 




2.319 xlO" 1 


2.318 xlO" 1 


1.101 xlO" 4 


3069 




s iMi-) 


2280 


1.869 xlO" 1 


1.863 xlO" 1 


5.418 xlO" 4 


2706 




" 3 P 2 "(2") 




2.100 xlO" 1 


2.099 xlO" 1 


1.203 xlO" 4 


2759 





Table 5. B{n — 2) meson mass spectra (first order). Units are in MeV. 



■ 2s+L Lj(J p ) 


Mo 


ci/Mo 


Pi /M 


ni /Mo 


Mcalc 


Mobs 


L So(0-) 


5849 


0.919 xlO' 2 


0.831 xlO" 2 


8.849 xlO" 4 


5902 




3 Si(i-) 




1.634 xlO" 2 


1.629 xlO" 2 


5.867 xlO" 5 


5944 




3 Po(0+) 


6025 


1.375 xlO" 2 


1.289 xlO" 2 


8.590 xlO" 4 


6108 




" 3 Pi"(l+) 




2.226 xlO" 2 


2.221 xlO" 2 


4.779 xlO" 5 


6160 




" 1 Pi"(l+) 


6098 


1.886 xlO" 2 


1.881 xlO" 2 


5.412 xlO" 5 


6213 




3 P 2 (2+) 




2.108 xlO" 2 


2.107 xlO" 2 


1.001 xlO" 5 


6227 




3 Dr(l-) 


5888 


1.610 xlO" 2 


1.605 xlO" 2 


4.668 xlO" 5 


5982 




" 3 P 2 "(2") 




1.809 xlO" 2 


1.808 xlO" 2 


1.037 xlO" 5 


5994 





Table 6. B s (n — 2) meson mass spectra (first order). Units are in MeV. 



' 2s+i Lj(J p ) 


Mi 


ci/M) 


Pi /M 


m/Mo 


Mcalc 


Mobs 


L S (Q-) 


5936 


0.878 xlO" 2 


0.802 xlO" 2 


7.588 xlO" 4 


5988 




3 Si(i-) 




1.597 xlO" 2 


1.592 xlO" 2 


5.013 xlO" 5 


6031 




3 Po(0+) 


6063 


1.399 xlO" 2 


1.325 xlO" 2 


7.429 xlO" 4 


6148 




" 3 Pi"(l+) 




2.299 xlO" 2 


2.292 xlO" 2 


6.908 xlO" 5 


6202 




" 1 Pi"(l+) 


6193 


1.828 xlO" 2 


1.826 xlO" 2 


2.351 xlO" 5 


6306 




3 P 2 (2+) 




2.052 xlO" 2 


2.052 xlO" 2 


4.709 xlO" 8 


6320 




3 Di(l-) 


5999 


7.631 xlO" 2 


7.627 xlO" 2 


4.449 xlO" 5 


6456 




" 3 P 2 "(2") 




-3.802 xlO" 2 


-3.802 xlO" 2 


4.861 xlO" 8 


5770 





